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We provide for the first time a complete parametrization for the matrix elements of the generic 
asymmetric, non-local and gauge-invariant canonical energy-momentum tensor, generalizing there¬ 
fore former works on the symmetric, local and gauge-invariant kinetic energy-momentum tensor 
also known as the Belinfante-Rosenfeld energy-momentum tensor. We discuss in detail the various 
constraints imposed by non-locality, linear and angular momentum conservation. We also derive the 
relations with two-parton generalized and transverse-momentum dependent distributions, clarifying 
what can be learned from the latter. In particular, we show explicitly that two-parton transverse- 
momentum dependent distributions cannot provide any model-independent information about the 
parton orbital angular momentum. On the way, we recover the Burkardt sum rule and obtain similar 
new sum rules for higher-twist distributions. 

PACS numbers: 11.15.-q,12.38. Aw,13.88.+e,13.60.Hb,14.20.Dh 


I. INTRODUCTION 

Following the canonical procedure based on Noether’s theorem, one ends up with a canonical energy-momentum 
tensor which is usually neither symmetric nor gauge invariant. Because of these pathologies, one often abandons 
the canonical energy-momentum tensor in favor of the Belinfante-Rosenfeld improved energy-momentum tensor |l]-0] 
which is both symmetric and gauge invariant. The Belinfante-Rosenfeld tensor differs from the canonical tensor by 
a so-called superpotential term which modifies the definition the momentum density but leaves both the total linear 
and angular momenta unchanged. It has the peculiar feature that it denies the mere existence of spin density. Indeed, 
from the conservation of total angular momentum d^J pvp = 0 where J p " p = r y T MP — r p T p " + S pvp with T pi/ the 
conserved total energy-momentum tensor and S p,vp the spin density tensor, one deduces that the antisymmetric part 
of the energy-momentum tensor is intimately related to the quark spin density T" p — T pv = —d ll S pMp . So, in the 
Belinfante-Rosenfeld approach, what is usually refered to as “spin” is simply described as a flow of momentum. There 
is therefore no clear distinction between spin and orbital angular momentum (OAM) in this approach, just like there 
is no clear distinction between energy flow T l ° and momentum density T° l . 

On the other hand, spin is an intrinsic property of a particle defined as one of the two Casimir invariants of the 
Poincare group (the other Casimir invariant being the mass). Contrary to OAM, one cannot change the spin of 
a particle by changing the Lorentz frame. Spin and OAM are distinguishable, and so dealing with a symmetric 
energy-momentum tensor is not very natural in Particle Physics. Where does this symmetry requirement come from? 
It is mainly motivated by General Relativity where gravity couples to a symmetric energy-momentum tensor. It is 
however important to notice that General Relativity is a classical theory while spin is fundamentally a quantum 
concept. Moreover, the symmetry of the energy-momentum tensor in General Relativity follows from the postulated 
absence of space-time torsion. More general theories relax the no-torsion assumption and do not require the energy- 
momentum tensor to be symmetric. The gravitational effects of the antisymmetric part of the energy-momentum 
tensor are however extremely small and are expected to show up only under extreme conditions, see e.g. [Ji6| and ref¬ 
erences therein. Finally, we note that the classical argument in favor of a symmetric energy-momentum tensor based 
on dimensional analysis and presented e.g. in section 5.7 of [7], is valid only for the orbital form of angular momentum. 

The early papers about the proton spin decomposition 0-03 start with the Belinfante-Rosenfeld tensor, but then 
add appropriate superpotential terms to decompose the quark angular momentum into spin and orbital contributions. 
According to textbooks EMI, no such decomposition is possible for the gauge field angular momentum. Though, 
photon spin and OAM are routinely measured in Quantum ElectroDynamics, see e.g. [131 ] and references therein. In 
Quantum ChromoDynamics (QCD), a quantity called AG which can be interpreted in the light-front gauge A°+A 3 = 0 
as the gluon spin [8|] has been measured in polarized deep inelastic and proton-proton scatterings, see fl4 ] for a recent 
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TABLE I: Properties of the various forms of the energy-momentum tensor in a gauge theory. 


Family 

Energy-momentum tensor 

Gauge invariant 

Local 

Symmetric 


Belinfante-Rosenfeld [G3] 

/ 

/ 

/ 

Kinetic 

Ji [9] 

/ 

/ 

- 


Wakamatsu [4L] 

/ 

- 

- 

Canonical 

Jaffe-Manohar [5| 

- 

/ 

- 


Chen et al. [15] 

/ 

- 

- 


analysis. In order to account for these experimental facts, Chen et al. claimed in 2008 that the textbooks were wrong, 
and proposed a formal gauge-invariant decomposition of the photon and gluon angular momentum 15]. This triggered 
a lot of criticism and an outpouring of theoretical papers, summarized in the recent reviews fliil. Il7l] . The apparent 
contradiction with the textbook claim was solved by realizing that the Chen et al. construction is intrinsically non¬ 
local Il8l420j ] , whereas textbooks implicitly refered to local quantities only. It has actually been known for quite some 
time that gauge invariance can be restored by allowing the quantities to be non-local [2l|, |22[ • Although there are in 
principle infinitely many gauge-invariant non-local quantities reducing formally to the same gauge non-invariant local 
expression in the appropriate gauge, the experimental conditions ultimately determine which ones are accessible [23[ . 

Parton distributions are typical examples of measurable non-local quantities. Gauge invariance is ensured by 
a Wilson line whose path is determine by the factorization theorems [24j. Ji has shown that the kinetic OAM, 
which is local and gauge invariant, can be expressed in terms of Generalized Parton Distributions (GPDs) that are 
accessible in some exclusive experiments like e.g. Deeply Virtual Compton Scattering [9]]. Since the local expression 
for the canonical OAM is gauge non-invariant [§}, it was thought for a long time that it cannot be measured 
and should therefore be considered as unphysical. The situation has changed once it has been realized that the 
non-local expression for the canonical OAM, which is gauge invariant, can be expressed in terms of fcr-dependent 
GPDs, also known as Generalized Transverse-Momentum dependent Distributions (GTMDs) [25l - l27j . GTMDs 
are extremely interesting as they provide the maximal information about the phase-space or Wigner distribution 
of quarks and gluons. Unfortunately, apart possibly in the low-a; regime, it is not known so far how to access 
these GTMDs experimentally [28]. The situation is however not hopeless since GTMDs can be accessed indirectly 
using realistic models, see e.a. [25l. l29U34| . Another possibility is to compute the GTMDs on the lattice. The 
traditional approach is to compute moments of the parton distributions using a tower of gauge-invariant local 
operators. Unfortunately, this approach does not allow one to compute AG because the latter does not correspond 
to any gauge-invariant local operator. However, new strategies have recently been proposed allowing in principle 
the computation of matrix element of non-local operators on the lattice, and show already encouraging results |35|-|40|. 

Many different forms have been proposed for the energy-momentum tensor in a g auge theory. Their properties 
are summarized in Table [T] All the forms can be sorted into two families [l9j, |4l], l42j : kinetic (or mechanical) and 
canonical. They all give the same total linear momentum, but attribute different momentum densities to the various 
constituents. To the best of our knowledge, only the matrix elements of the local energy-momentum tensors have 
been discussed in the literature so far. The first complete parametrization of the matrix elements of the symmetric 
gauge-invariant local operator (i.e. Belinfante-Rosenfeld tensor) has been given in |j§] and further discussed in J43|. 
The matrix elements of the asymmetric local gauge-invariant operator (i.e. Ji tensor) has first been discussed 
in 0, but the correct parametrization in the off-forward case was given in [44| . There has also been a simple 
attempt to parametrize in a similar way the matrix elements of the asymmetric local gauge non-invariant operator 
(i.e. Jaffe-Manohar tensor), but this led to the absurd conclusion that canonical and kinetic matrix elements are 
the same [45J. We will argue in the present paper that the failure of this attempt can be understood as due to 
the absence of an important piece of information in the parametrization of [45j . We will also show explicitly that 
two-parton Transverse-Momentum Distributions (TMDs), though sensitive to OAM, cannot provide any quantitative 
model-independent information about the OAM. 

The paper is organized as follows. In sectionQj] we decompose the QCD energy-momentum and generalized angular 
momentum tensors into quark and gluon contributions, and compare the various forms found in the literature. In 
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section Hill we provide for the first time the parametrization of the generic non-local light-front gauge-invariant energy- 
momentum tensor and discuss various constraints in section llVl In section[Vl we derive the relations between the scalar 
functions appearing in this parametrization and derive the relations with the two-parton generalized and transverse- 
momentum dependent distributions, obtaining on the way new sum rules. Finally, we gather our conclusions in 
section IVII Some details about the parametrization are given in Appendix [A] 


II. THE GAUGE-INVARIANT LINEAR AND ANGULAR MOMENTUM TENSORS 


In order to deal most conveniently with the various gauge-invariant decompositions proposed in the literature, we 
consider the following five gauge-invariant energy-momentum tensors 

Tf"(r) = i>(r)'y' i %D v ip(r), 

rrW = —2Tr [G pa (r)G\(r)) + g pv ±Tr[G^(r)G^(r)], 

T r ( r ) = (1) 

Tf» = \e^d a \${r) m ifr(r)], 

Tr(r) = -2d a Tr[G» a (r)A" hy s(r)\, 

AA . AA AA —y -4 — 

where £0123 = +1 and + gAf 1 is the hermitian covariant derivative with 8^ = 8^ — 8^. Similarly, we 

consider the following seven gauge-invariant generalized angular momentum tensors 


K vp {r) 

S pvp (r) 

sr(r) 


= r v T pp (r) — r p T pv (r), 
= i e ^ CT ^(r)7 <T 7 5 V>(r), 
= -2Tr[G'*[*'(r)^L(r)' 


a = !,••• ,5, 


( 2 ) 


where x^y 1 ^ = x p y v — x v y p . The standard expressions for the Belinfante-Rosenfeld, Ji, Wakamatsu and Chen et al. 
decompositions 1 are then obtained by combining these contributions according to Tables |TT] and Mil and using the 
following identities based on the QCD equations of motion 


= -e pua,3 d a \4>(r)'yp'y 5 i’{r)\ , ^ 

2[D a G af) (r)] C c , = -g^(r) 7 ^f(r), 

where c, d are color indices in the fundamental representation and = 8^ — ig[A^, } is the adjoint covariant 
derivative. In particular, because of the first identity in Eq. (Jd)), we can write T pv (r) = — \ (r) and therefore 

discard the tensor T pv (r) in the following discussions. Note that the tensors T\ pu \r), T pu (r) 1 L plJp {r) — S pvp (r) and 
L pvp (r) + S 2 Vp {r) have the form of a superpotential 8 a f^ ap ^'"(r) Assuming as usual that surface terms vanish, 
this means that we have 


d^ v] (r) = 0, 

J d 3 rrj H (r) =0, 

W(r) = 0, 

J d 3 r T™(r) = 0, 

lf pl (r) = -8^Sr P (r), 

J d 3 rL™ p {r)= J d 3 rSr P (r), 

Tt P] (r) = -d^Srir), 

J d 3 r L™ p {r) = — j d 3 rST P {r) 


(4) 


1 We used the original covariant form of Ref. [42 1 and not the one in Ref. El which differs only by how one separates the pure-boost 
terms into quark and gluon contributions. 
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TABLE II: Expressions for the Belinfante-Rosenfeld, Ji, Wakamatsu and Chen et al. forms energy momentum tensors for quarks and 
gluons. 



Belinfante-Rosenfeld 

Ji 

Wakamatsu (gik) 

Chen et al. (gic) 

Tq W {r) 

T^{r) + T pv {r) 

TTW 

T-Hr) 

T^{r) + Tg v (r) 

T^(r) 

ITM 

7TM 

TW + TW 

rrM-TTM + TTM 


TABLE III: Expressions for the Belinfante-Rosenfeld, Ji, Wakamatsu and Chen et al. forms of the generalized spin and orbital angular 
momentum tensors for quarks and gluons. 



Belinfante-Rosenfeld 

Ji 

Wakamatsu (gik) 

Chen et al. (gic) 


0 


Sf" p (r) 

ST P (r) 

L>f p (r) 

Lr p (r) + L^(r) 

LZ up {r) 

L^ p {r) 

L pl,p {r) + L p " p {r) 

srM 

0 

0 

sr(r) 

sr(r) 

Lg"(r) 

UT{r) 

Lf P W 

Lr(r)+Lr(r ) 

Lr(r)-Lr(r) + Lr(r) 


where n is a timelike or light like four-vector and d 3 r = € a p^s dr^ A dr 7 A d r s is the volume element. This ensures 
that the quark and gluon linear and angular momenta are the same in the three kinetic decompositions 



where J pvp (r) = S pvp {r) + L pvp (r). 


The Wakamatsu and Chen et al. decompositions require the introduction of a pure-gauge field 

J 4Pure( r ) = | W ( r )^ W - 1 ( r ) ) (6) 

where W(r) (called C/ pure (r) in 0) is some phase factor transforming as W(r) U(r)W(r) under gauge transfor¬ 
mations. The “physical” gluon field is then defined as 

Afy*(r) = A„(r) - A^r). ( 7 ) 

In the gauge where W(r) = 1, the Chen et al. decomposition takes the same mathematical form as the Jaffe-Manohar 
decomposition, and can therefore be considered as a gauge-invariant extension of the latter fl6[ fl9l l46l l47j . The 
phase factor W(r) is non-locally related to the field strength and is in principle not unique [161 Il9|. The original 
Wakamatsu [4lJ and Chen et al. '0 decompositions correspond to a particular choice of the phase factor which makes 
the physical field transverse in a given Lorentz frame. Leaving the phase factor unspecified allows us to consider at 
once two whole classes of decompositions differing simply by the precise form of the non-local phase factor. In order 
to stress this point, we will follow from now on the terminology of Ref. [l6| and refer to the Wakamatsu and Chen et 
al. decompositions as the gauge-invariant kinetic (gik) and canonical (gic) decompositions, respectively. For a given 
phase factor, the difference between the gauge-invariant kinetic and canonical decompositions lies in the separation 
of total linear and orbital angular momentum into quark and gluon contributions. This difference corresponds to 


M^ p (r) 




which are called potential linear and angular momentum tensors [41l . |42| , respectively. 


( 8 ) 
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III. PARAMETRIZATION 


In practice, since we want to relate the matrix elements of the gauge-invariant energy-momentum tensor to measur¬ 
able parton distributions, we choose the non-local phase factor W(r) to be a Wilson line W n (r, r 0 ) connecting a fixed 
reference point ro (usually taken at infinity) to the point of interest r. According to the factorization theorems [24| , 
these Wilson lines run essentially in a straight line along the light-front (LF) direction given by a lightlike four-vector 
n to the intermediate point r n = r ±oon, and then in the transverse direction to r 0 . In some sense, these Wilson lines 
can be viewed as a background gluon field generated by the hard part of the scattering. The Wilson line associated 
with the first part of the path 


Wn(r,r n )=V 


D —ig f±°° n-A(r+Xn) dA 


(9) 


is responsible for making the LF gauge n ■ A = 0 special, since this is the gauge where W„(r, r n ) = 1. The transverse 
Wilson line Wn( r n, ro) is associated with the residual gauge freedom and can be set to 1 using appropriate boundary 
conditions for the gauge field [20l [27j ■ Our gauge-invariant canonical energy-momentum tensor will then be physically 
equivalent to the Jaffe-Manohar tensor considered in the LF gauge n • A = 0 with appropriate boundary conditions. 

We will consider in the following the generic LF gauge-invariant energy-momentum tensor of which the Belinfante- 
Rosenfeld, Ji, gauge-invariant kinetic and canonical energy-momentum tensors represent particular cases. The matrix 
elements of the generic LF gauge-invariant energy-momentum tensor depends in principle on n. This dependence 
was overlooked in |45j|, leading to absurd conclusions. Since any rescaled lightlike four-vector an specifies the same 
LF Wilson line, the matrix elements of the generic LF gauge-invariant energy-momentum tensor actually depends, 
beside the average target momentum P = (p 1 +p)/2 and the momentum transfer A = p' — p 7 also on the following 
four-vector 


N = 


M 2 n 
P ■ n 


( 10 ) 


with M the target mass, and on the parameter 77 = ±1 indicating whether the LF Wilson lines are future-pointing 
( r) = + 1 ) or past-pointing (77 = — 1 ). Note that, contrary to n , the lightlike four-vector N has the same dimension 
and transformation properties under space-time symmetries as the momentum variables. Since P ■ A = 0 and 
M 2 = P ■ N = P 2 + A 2 /4, the scalar functions parametrizing the matrix elements of the generic LF gauge-invariant 
energy-momentum tensor are functions of the two scalar variables £ = —(A • N)/2(P ■ N ) and t = A 2 . Choosing 
the standard form for the lightlike four-vector n = (1,0, 0,-1) leads to the usual expression £ = —A + /2 P + with 
a ± = a° ± a 3 . Because these scalar functions also depend on the parameter 77 , they are complex-valued just like the 
GTMDs @,0. 


Using the techniques from the Appendix A of Ref. [28|, we find that the matrix elements of the generic LF gauge- 
invariant energy-momentum tensor for a spin- 1/2 target can be parametrized as 

(p',S'lTr( 0 )\p,S) = u(p',S')r^(P,A,N; V )u(p,S), ( 11 ) 

where S and S' are the initial and final target polarization four-vectors satisfying p-S = p'-S 1 = 0 and S 2 = S' 2 = — M 2 , 
and T ^ stands for 


= Mg^AI + 
N»N V 

Mg^ v Bl - 
P^A 


pH pv 

_ A a 4 

M ^ 24 
P /7 

Bn 


A M A" 




P^ia uA 


M ‘ 3 




+ 


+ 


+ 


M 
P^ia vN 
2 M 
P^A" 


£*14 


M 

ptipv 

M 

P"A^ 


B a a 


M 

A^A" 


B-i + 


2 M ' 4 
N^ia vA 


Al 


2 M 


B% 


B. 


oa 

M 

p-'i^N 


M 10 
A 11 N v 


N^N 1 ' „ 
B“ + B 


M 


pa 

^16 


M 11 
A 


P v ia^ A 
2 M 

N v itj^ A 
2 M 

piiKTV 

B 


■B? + 


A »io vN 


2 M 


■B, 


A v ia tiN 


p u A - / 1 


M 


12 


Dfl 

B 13 


M 


B 


17 


M 

— i<T» v B“ 


18 ' 


M 
A v ia^ A 
' 2 M 


2 M 


B a 7 


N A 


2M 2 


( 12 ) 


Dfl 

-°19 


20 


M 


■B: 


2 M 
P U A^ 


B° 


N^ia uN 


2 M 


Bn 


N v ia^ N 


2 M 


II 


24 


II 


a (, r 


M 25 ' M M ' M 

For convenience, we used the notation ia^ b = ia^ a b a . The factors of i have been chosen such that the real part of 
the scalar functions is 77 -even and the imaginary part is 77 -odd 

A“(£,t;77) = X e ’ a (£,t) + 777X ? °’ a (£,t) 


II 


A^AT/7 


B. 


27 


1(7 


23 
N A 


2 M 2 ' 


(13) 
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as a consequence of naive time-reversal symmetry. Hermiticity then implies that the real part of Bj with j > 14 is 
£-odd and the imaginary part is £-even. For the other functions, the real part is £-even and the imaginary part is 
£-odd. 

We have found that the parametrization of the matrix elements of the generic LF gauge-invariant energy-momentum 
tensor for a spin-1/2 target involves 32 complex-valued scalar functions. This number can be obtained from a simple 
counting. The generic energy-momentum tensor T'g v has 4 x 4 = 16 components. The target state polarizations ±S 
and ±S' bring another factor of 2 x 2 = 4, but parity symmetry reduces the number of independent polarization 
configurations by a factor 2, leading to a total of 32 independent complex-valued amplitudes {p', S'\T^ (0)\p, S). These 
32 independent amplitudes correspond to 32 independent Dirac structures, a particular set being given by Eq. m- 
Any other Dirac structure like e.g. ia^ p or i e ^ NA ^ can be expressed onshell as a linear combination of these 
32 structures, see Appendix lAl of this paper. 


IV. CONSTRAINTS 

The parametrization (1121) is very general and does not take into account several constraints like linear and angular 
momentum conservation. We discuss in this section the various constraints and the relation to former works on the 
local gauge-invariant energy-momentum tensor. 

For latter convenience, we introduce the Sudakov decomposition of a generic four-vector 


= (a • n)n M + (a • n)n M + a(/. 


(14) 


together with the transverse Kronecker and Levi-Civita symbols 


6£ v = SjJ - n^n v - 


ett!' = e^ al3 7 


'0^/3 1 


(15) 


where n is the lightlike four-vector satisfying n ■ n = 1 and such that = 0 . 


A. Local operators 


The energy-momentum tensors T pv (r) and T 2 V (r) are local. The corresponding matrix elements cannot therefore 
depend on N or g. All the scalar functions must then vanish except the five real-valued functions A £ j’ a (0,t) with 
a = 1,2. These are related to the standard (local) energy-momentum form factors (FFs) £1, HE El as follows 


A q {t) = A e 2 A (0,t), 

Adt) = A 2 2 (0,t ), 


B q (t) = Af(0,t) + Af(0,t) - Af(0,t), 

Bait) = Af( 0, t) + Af (0, t) - Af( 0, t), 


C q (t) = Af(0,t), 

C G (t) = Af(0,t), 

(16) 

C q (t) = Aj’ x ( 0, t) + -Ay Af(0, t), 

C G (t) = Af(0,t) + 1 ^Af(0,t), 


D q (t) = Af (OA)-Af (0,f), 

0 = Af(0,t)-Af(0,t). 


The first four form factors parametrize the symmetric part of the local gauge-invariant energy-momentum tensor, 
whereas the last one parametrizes its antisymmetric part. Since T 2 v (r) is symmetric, we have Al’ 2 (0,t) = ^s’ 2 (0,t)- 


From our choice of the phase factor 


B. Light-front constraints 

it follows that A p h ys ■ N = o [3 [S3, leading to 


T 3 N ( r ) = T£ N (r) = 0 . 


(17) 
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Contracting our generic parametrization (1121) with N v , we find the relations 

A a 2 (U)-2£BU(Z,t)=0, 

-2ZA a 3 (ZA) + B? 5 (ZA)=0, 
t ) + B a (Z, t) - 2ZB a 24 (Z, t) = 0 , 

A a 5 (ZA)~^ZBUZA) = 0, 

BUZ, t) + BUZ, t ) + B“ 3 (Z, t) - 2 ZB^iZ, t) = 0, 

-2 ZBUZ, t ) + Bf 8 {Z, t) + BU(Z, t) = 0, 

-2ZBU(Z,t) + BU(Z,t)=0, 

for a = 3, 5 which we refer to as the LF constraints. 

C. Four-momentum conservation 

The total energy-momentum tensor T M ^(r) = T{ t! ' (r) + T 2 v (r) and the superpotential terms T\^ v \r) and T 3 v (r) 
are all conserved d fl T liV (r) = = d ll T^ L '(r) = 0. This translates at the level of the matrix elements as 

A^{p',S'\T^mP,S)=0, 

AUp , ,S , \T l 1 ^ ] (0)\p,S)=0, (19) 

AUp',S'\Tr(0)\p,S)=0, 

and implies the following constraints 

j2[Arm+^Arm] = ^ aw = 

a— 1 a—q,G 

A 5 1 (Z,t) + 1 fcAUZ,t)-2ZB 5 17 (Z,t) = 

-2ZB\{ZA) + BUZA) = 

-2ZBUZA) + ^ BUZA) = 

- 2 ZBKZA) ~ B 5 18 (ZA) = 

T^Bl(ZA)-2ZB 5 22 (ZA) = 

—b 7 (Z , t) + Bg(Z,t) + jjp b 30 {Z,t) - ^ZB 27 (Z,t ) = 

—2ZB 31 (Z, t) ~ B 23 (Z > t) + JT? B 2 q(Z11) = 

-2ZB 5 13 (Z, t) - B 5 21 (Z, t) + AL, B 25 (Z, t ) = 

which are compatible with Eq. m- 


0 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 , 


( 20 ) 


D. Forward limit and momentum 


In the forward limit A —> 0, the parametrization of the generic LF gauge-invariant energy-momentum tensor reduces 


to 


(P, 5 , |T q ^(0)|P, S) = 2MV'X’ Q + 2P^P u A 2 ’ a + 2 N^N u B 3 ' a + 2P^N v B e 2 a + 2P , 'N^B e 3 a 


+ t) e^ SP B°’ s a - 


V 


vS ~qO,cl 
20 


P^e v T s B. 


■ P^etf B°[ a 


+ B 22 a + N v dj? B. 


1/ ^_flS -j^O)CL 

23 


( 21 ) 


Since T§ v (r) is a total divergence, its matrix elements are proportional to A and therefore vanish in the forward limit, 
leading to 


Af(0,0) = Af(0,0) = 0, 

P 4 ’ 5 ( 0 , 0 ) = P 2 ’ 5 ( 0 , 0 ) = Pg’ 5 ( 0 , 0 ) = 0 , 


(22) 


Big (0- 0) = P 2 ° o °(0,0) = P 2 °i°(0,0) = P 22 °(0,0) = P 2 °3°(0,0) = 0. 


,o,5/ 


^0,5 


o,5 / 


,o,5/ 
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Moreover, since the tensors T( 11/ (r) and T 2 u (r) are local, the only non-vanishing scalars £?“(0,0) arise from the 
potential term T pu {r). This means in particular that naive T-odd effects in the forward limit are necessarily 
associated with the canonical momentum and disappear when summed over all partons. 


Contracting now Eq. (12T1) with A), gives the average four-momentum in the LF form of dynamics 

<P«> = ^ (P, S\T? V (P)\P, S) = P»Ar + N"(Ar + B e 2 ’ a ) + f ef (B& - B%). 


(23) 


In particular, using Eq. (1161) we recover the standard expression for the gauge-invariant kinetic four-momentum in 
terms of the energy-momentum FFs 


CPgik.J = P u A a (0) + N v C a { 0), a = q,G. (24) 

Interestingly, the last term in Eq. (l23| is naive T-odd and can be interpreted as the spin-dependent contribution to 
the momentum arising from initial and/or final-state interactions. Because of the structure this naive T-odd 
contribution is transverse and requires a transverse target polarization. As we will see in section fVBl this is related to 
the Sivers effect [49]. The combination of scalars A®’“(0, 0) + B 2 ’ a '(Q, 0) contributes only to the energy and is therefore 
related to the interaction term in the Hamiltonian. In the forward limit, the LF constraints (1181) imply that 

Af(0,0) = 0, (25) 

unless Bif(£,t) behaves as l/£ near £ = 0. This suggests that the scalars A 2 ( 0,0) and A2 ,2 (0,0), and hence A g (0) 
and Ag(0 ), can be interpreted as the interaction-independent contributions of, respectively, quarks and gluons to 
the four-momentum. This is further supported by the observation that ^’“(O,!)) is the only contribution to the 
longitudinal momentum (p™), which is purely kinematical in LF quantization. In other words, there is no difference 
between the longitudinal component of the average kinetic and canonical momenta. 

Finally, since the total four-momentum is (p v ) = P v , we obtain from Eq. (1231) the momentum constraints 


Y ^’ 1 (o,o)= Y £*(°) = o- 

i=l,2 a=q,G 

Y A e /(0,0)= Y 4,(0) = !, 


(26) 


which are consistent with Eq. (1201) . In particular, the vanishing of the average total transverse momentum, known 
as the Burkardt sum rule [50|, |51[ , is trivially taken into account in our parametrization because the potential term 
To*' 1 ' (r), and hence B/g 3 (0, 0) — I?2o 3 (0,0), drops out of the sum over all partons. 


E. Angular momentum 

Since we have a complete parametrization of the matrix elements of the generic LF gauge-invariant energy- 
momentum tensor, we can easily compute the matrix elements of the corresponding OAM tensor L pup (r) given 
by Eq. ©. Because of the explicit factors of position r, the matrix elements of the generic LF gauge-invariant OAM 
tensor need to be handled with care [3,01. Focusing on the longitudinal component of OAM, we find 

. ULX CX J A=0 

= j^r A 4 “(0,0). (27) 

For a longitudinally polarized target, we have S ■ N = M 2 and so A%’ a (0, 0) can be interpreted as the average fraction 
of target longitudinal angular momentum carried by the OAM associated with the energy-momentum tensor T pu (r) 
in the LF form of d yna mics. We confirm in particular that the integrated OAM does not receive any naive T-odd 
contribution [20, [23, Ell [S3 ■ Using Eq. (fl6l) . we also recover the standard expressions for the Bclinfante and Ji forms 
of longitudinal OAM in terms of the energy-momentum FFs [H, [l(J El 04] 


(LI) = ^ {P, S\\e TaP L N a ^mP, S) = ^ 2 


(^Bel.i) = <£&el,L> = 5 [4,(0) + £ o (0)] , 

(L^l) = | [4(0) + B q (0) + D q (0)} ^ 


a = q,G, 


(28) 
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Remarkably, thanks to Eq. 0 we can also express the quark and gluon spin contributions in terms of the scalar 
functions parametrizing the generic LF gauge-invariant energy-momentum tensor. From the integral relations in 
Eq. 0, we find that the quark and gluon longitudinal spin contributions are given by 


(51) ee (P,5||e T ^5f^(0)|P,5) 
(S?) = ^2 (P,S\^TafiS^(0)\P,S) 



Af(0,0)-Af(0,0) 


S-N 

AP 


Af ( 0 , 0 ). 


S-N 

Up-’ 


(29) 


where we have used Eq. 0 to express L pvp {r) in terms of T^ v (r). The scalars —\\A!± 1 ( 0,0) — Ag’^O, 0)] = — \D q (d) 
and —A e p‘(Q 1 0) can therefore be interpreted as the average fraction of target longitudinal angular momentum carried 
by the spin of quarks and gluons, respectively. It is easy to check that the following relations 


(S a L ) + (L a g ^ L ) = (J^ hL ), a = q,G (30) 

are satisfied and that the longitudinal component of the potential OAM is given by 

q ]\f 

(L Q eiC , L ) ~ (L^ >l ) = - [{L°. l ) - (L% >l )] = Af( 0,0). (31) 

Note that the differential relations in Eq. 0, which translate at the level of matrix elements as 

(p',5'|T'^ ] (0)b,5) = -iA,{p\S'\SrmP,S), 

(p',5'|Tf p] (0)|p,5) = -*A M (p',5'|5r P (0)b,5), 

do not provide additional constraints. Indeed, at O(A 0 ), they just reduce to the antisymmetric part of the forward 
limit (ED- At higher orders in A, the identification of coefficients between the LHS and the RHS of Eq. (1321) are 
spoiled by the condition P • A = 0 which follows from the onshcll relation for the target (P ± y) 2 = M 2 . 

Finally, since the total angular momentum is 1/2, we obtain from Eqs. (1271) and (1291) the angular momentum 
constraints 


E K“(0,0) + Al’ a (0 ,0)] = E [M0) + a»(0)] = 1. (33) 

a—1,2 a—q,G 

Combined with the momentum constraints (1261) . this leads to 

E [At’ a (0, 0) + AJ Q (0,0) - A e 2 ’ a (0, 0)] = E S -(°) = (34) 

a=l,2 a—q,G 

known as the anomalous gravitomagnetic moment sum rule l53l. [54] . 


V. LINK WITH MEASURABLE PARTON DISTRIBUTIONS 

Now we are going to see how the scalar functions parametrizing the matrix elements of the generic LF gauge- 
invariant energy-momentum tensor are related to GPDs accessed in exclusive scatterings [55j and TMDs accessed in 
semi-inclusive scatterings J24|. For convenience, we shall focus in the following on the quark sector. The gluon sector 
proceeds analogously. 


A. Generalized Parton Distributions 

The quark vector GPD correlator is defined as 

F S'S (P,U A, N) = (P -n) J ^ e ix A p - n ) ( p ', ±n)^n)\p, S). (35) 

Remarkably, its second Mellin moment is related to the matrix elements of the quark LF gauge-invariant energy- 
momentum tensor [Til Hi], |55| 

J dxx FpSJ (P, x, A, N) = ^ (p' ) S'|T'‘"(0)|p,S>. 


(36) 
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Note that we do need to specify whether this corresponds to the kinetic or canonical version of the LF gauge-invariant 
energy-momentum tensor simply because T^ q (r) = " (r) owing to Eq. (1171) . 

Up to twist 4, the quark vector GPD correlator (l35l) is parametrized as [H] 


FW{P,x,A,N) 

F [ ;,t ] (P, X ,A,N) 

FP s (P,x,A,N) 


= u(p',S') 

= u(p',S') 

M 2 

“ (P • n) 2 


jiH q 


-71A 


2 M 


■E q 


u(p, S), 

fiA “ 


MialF 1 „ 


P ■ n 
u(p\ S' 


H q 2T - Yt (e! t ~ dE q T ) 


i(p, S ), 


mi 


io 


nA 


E q 

2 M 3 


u(p,S). 


(37) 


From Eq. we then find the relations between the second Mellin moment of vector GPDs and the energy- 

momentum FFs in the quark sector 


J dxx H q (x,Z,t) 
j d xx E q (x,Z,t) 
J d x x H q T (x,Z,t) 

/ d xxE 2T (x,Z,t) 


— Aq(t) + 4£ 2 Cq(t), 

= B q {t)-mC q {t), 
= 0 , 

= 0 , 


(38) 


J d xxH q T {x,Z,t) = -2 £C q (t), 

J dxxE q T (x,Z,t) = —\[A q (t) + B q {t) — D q (t)\, 

J d XX H$(x,£,t) = \A q {t)+C q {t) - 2 £ 2 J^Cq{t) + -Ap- [Bq(t) - 8 Cq(t) - Dq(t)\, 

J d xx[H%{x,Z,t) + E${x,Z,t)] = 2 

The relations involving twist-3 GPDs are consistent with those found in Ref. [56] where the parametrization is related 
to the one we used as follows 


H% T {x, Z, t ) = —2£ G|(x, Z, t), 

E q 2T (x, £,t) = 2 [G q (x, Z, 1) - £G\{x, £, *)], 

H q 2T {x. Z, t ) = \G\(x, Z, t), 39 

H q (x, Z, t) + E q (x, Z, t) + E q T (x, Z, t) = -G q (x, Z, t ) + 2[ZG q 3 (x, Z, t) - G q (x, Z, *)]• 

The explicit relations involving twist-4 GPDs are new but somewhat academical as these functions are much harder 
to access experimentally 2 . 

Thanks to Eq. (1321) . the antisymmetric part of the LF gauge-invariant kinetic energy-momentum tensor can be 
related to the local axial-vector correlator 

- |e^ AQ | dxF^(P,x,A,N) = {p\S'\T [ g ^ q m P ,S). (40) 

From the parametrization f55j ] 


d xF [ f s ^ ] 


(P,x,A,N) = u(p',S') 


7 M 75 G^(f) 


A^7s 
2 M 


Gp(t) 


i(p, S), 


(41) 


2 We note in passing a typo in Ref. M where a factor ^ is missing in front of the D q (t) energy-momentum FF in the RHS of Eq. (432). 
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where G q A [t) = f dx H q (x, £, t) is the axial-vector FF and G q P (t) = f dx E q (x,t;,t) is the induced pseudoscalar FF, it 
is easy to show using onshell identities that [13, [IE, 0 

G\{t) = —D q (t). (42) 

This is naturally consistent with the observation in section [IVE] that the scalar ~^D q ( 0) can be regarded as the quark 
spin contribution to the total angular momentum in a longitudinally polarized target. 


B. Transverse-Momentum dependent Distributions 

The quark vector TMD correlator is defined as 3 

^(P,x,k T ,N-rj) = (P-n) J d (fc-n) J ^e ifc -*<p ) S , |?(-f)7' 1 W n (-f,fMf)|p,S). (43) 

Similarly to the GPD case, the second Mellin moment of this correlator is related to the forward matrix elements of 
the quark LF gauge-invariant energy-momentum tensor 

J dxd 2 k T x^(P,x,k T ,N-, V ) = ^ {p,S'\T£ N (0)\p,S). (44) 

Remarkably, treating with due care the LF Wilson line 0, HE, 0 SS|, one can similarly show that the second 
transverse moment of J is related to the forward matrix elements of the quark LF gauge-invariant canonical 
energy-momentum tensor 


J dx d 2 kr k P $g,J(P, x, kr, N; rj) 


^ v {p,s'\Tr q (o)\ P ,s). 


(45) 


That the relation holds for the canonical version of the LF gauge-invariant energy-momentum tensor is determined 
by the particular shape © of the Wilson line. Working instead with a straight Wilson line connecting directly 
the points the relation (1451) would hold for the kinetic version of the LF gauge-invariant energy-momentum 
tensor 0, [63, (63 • We stress once again that our choice for the Wilson line © was simply motivated by the fact that 
factorization theorems require Wilson lines that run essentially along the LF direction n [IE] . 

Up to twist 4, the quark vector TMD correlator (H3l) is parametrized as [H, [63:, 0 


*fk p > k T ,N;r,) = 2(P- n) f\ - r, fa f# 


P,x,k T ,N-, V ) = 2M 


S. f-L<7 ^ e T fQ _ 

M J '/ M JT 


nS 


d>f ] s (P,x,k T ,N;ri) = 


P ■ n i 


M j 'I M -IT V 

f: 


{k^kTv n e T^ ( S-n ) 

---/t - V -—- 






(P-n)M JL 


ff 1 


(46) 


2M ' r /l - V ^ f±q 


M 1 J 3T 


where we have extracted explicitly the rj- dependence. From Eq. (TCT) and the forward limit A —> 0 of the LF 
constraints m, we find the relations between the second Mellin moment of vector TMDs and the energy-momentum 
FFs in the quark sector 


Jdxd 2 k T xmx,k 2 T ) = 4j(0), 


J dx d 2 kr x f P (x,k P ) = 0, 

J dxd 2 k T x f q (x, k P ) = \A q (d) + C q (0). 


(47) 


3 For simplicity, we considered the naive definition of the TMD correlator where the soft factor is not included [24|. This allows one to 
treat in a simple way the fc^-integrations. A more careful treatment based on the proper definition of TMDs proceeds analogously. 
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Together with Eq. (l38l) . these relations are consistent with the fact that the GPD and TMD correlators have the same 
collinear forward limit 


/ d 2 fcr (P, x, k T , TV; 77 ) = F l ^ ] (P, x,A = 0, TV) 


(48) 


which implies 


J d 2 k T f?(x,k$) = Hi(x,0,0), 

j d 2 fcr fcfO = 0 = iJ| T (a:, 0,0), (49) 

J d 2 k T fi(x,k*) = Hi(x,0,0). 


More interesting are the relations involving the second transverse moment of TMDs. From Eq. fl5l) and the forward 
limit A —> 0 of the LF constraints dlST) . we find 


dz d 2 fc T /#( x, k 2 ) = -I [p°g 3 (o, 0) - B°’ 0 3 ( 0,0) 


2M 

k: 


dx cl 2 k T f ±ll (x, k z T ) = C q { 0) + A e *( 0,0) 


e,3 / 


J dxd 2 k T 

J da; d 2 k T (x, k%) 


iBif(0,0), 

3 kf(0,0) + p 2 °6 3 ( 0,0) + 2B 2 °’ 3 (0,0)' 


(50) 


Anticipating that similar results hold in the gluon sector, we sum over all partons and obtain the following sum rules 



(51) 


The first sum rule is known as the Burkardt sum rule |50l I 5 H and simply expresses the fact that the total momentum 
transverse to the target momentum has to vanish. The other three sum rules are to the best of our knowledge new. 
They express the fact that the total flow of transverse momentum has also to vanish. They involve higher-twist TMDs 
and are therefore much harder to test experimentally. Nevertheless, it would be very interesting to test them using 
phenomenological models, Lattice QCD and perturbative QCD. 

As a final remark, we would like to stress that the above results show explicitly that TMDs cannot provide any 
information about the scalar ^’“(0, 0), which means no quantitative model-independent information about the parton 
OAM, as anticipated e.g. in j65| 


VI. CONCLUSIONS 

There has been a prejudice against the canonical form of the quark and gluon energy-momentum tensor, and 
consequently of the corresponding linear and orbital angular momenta, due to the fact that it cannot be written 
locally in a gauge-invariant way. A gauge-invariant expression can however be obtained by relaxing the locality 
requirement in a way that does not harm causality. This indicates that the canonical energy-momentum tensor can 
be considered as a physical object and measured experimentally. In particular, it can be accessed via particular 
moments of two- and three-parton correlators which are extracted from numerous physical processes. 
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In this study, we provided for the first time a complete parametrization for the matrix elements of the generic 
asymmetric, non-local and gauge-invariant canonical energy-momentum tensor. We found that a generic canonical 
energy-momentum tensor for a spin-1/2 target consists in 32 independent complex amplitudes. We discussed in detail 
the various constraints on these amplitudes imposed by non-locality, linear and angular momentum conservation. 
This generalizes therefore former works on the symmetric, local and gauge-invariant kinetic energy-momentum tensor 
also known as the Belinfante-Rosenfeld energy-momentum tensor. 

We also showed that some of the amplitudes can be expressed in terms of particular moments of two-parton 
generalized and transverse-momentum dependent distributions, and are therefore clearly measurable. In particular, 
we proved explicitly that two-parton transverse-momentum dependent distributions cannot provide any quantitative 
model-independent information about the parton orbital angular momentum. On the way, we recovered the Burkardt 
sum rule, expressing basically conservation of transverse momentum, and derived three new sum rules invvolving 
higher-twist distributions. We obtained these results by choosing the non-local phase factors defined by a lightlike 
four-vector n, in order to make contact with parton physics and factorization theorems. 

We believe the present paper will help clarify the differences between canonical and kinetic energy-momentum 
tensors, and their links with parton distributions. We also expect getting more insights into these matters in a near 
future coming from explicit results obtained within covariant models, Lattice QCD and perturbative QCD. 
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Appendix A: Parametrization 

From the discrete space-time symmetries, we find that the Dirac structure T ^ associated with the matrix elements 
of the generic LF gauge-invariant energy-momentum tensor has to satisfy the following constraints 


rre A > v) = 7 °r^(P, -A, TV; 77 ) 7 ° Hermiticity 

= 7°rg p (P, A, TV; 77)7° Parity 

= (— i"f 5 C)T* p/ '' (P , A, TV; — ?y)(—T 75 C) Time-reversal 

where C is the charge conjugation matrix and b = TP = ( 6 °, —b). Using the Gordon identities 


(Al) 


u{p\ S’)^u(p, S) = u{p\ S') 
0 = u{p. S’) 


PM i a v A 

+ 


M 2M 
A M ia^ p 


2 M M 


u{p , S), 
u{p, S), 


(A2) 


we can discard the 7 ^ and icr^ p structures from the parametrization. Similarly, we can discard the structure e lluaP ^5 
thanks to the following onshell identity 


u(p', S') [2e^ QP 75 + N Jj a va + A v o a » + A a a pv ] u(p, S) = 0, 

Contracting the e identity 

gaP^vpa ^appipc r/3 ^av^app. _|_ g°tp e >_|_ gUa^pvp _ q 

with PpApNa gives 

pa^pvAN ga[p e u]ANP _|_ ^a^pvNP _ jya^pvAP _ q 

Contracting further with N a: P a and A a leads to 

(P • N) e pvAN + N^e^ ANP + (A • TV) e pvNP = 0, 

p2 e pvAN _|_ p[p e u]ANP — (p . TV) e P vAP _ 

^[p e v]ANP + ^2 e pvNP _ (A • TV) e P vAP = 0. 


(A3) 

(A4) 

(A5) 

(A6) 
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Multiplying now by u{p', S')^u(j>, S) and using Eq. (IA3I) . the first two identities allow us to discard the structures 
e^ AN "f 5 and A^a v ^ A while the last identity is trivially satisfied. We are then left with the 32 independent structures 
given in Eq. ra¬ 
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